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Abstract 

In restored peatlands, recovery of carbon assimilation by peat-forming plants is 
a prerequisite for the recovery of ecosystem functioning. Restoration by rewett- 
ing may affect moss photosynthesis and respiration directly and/or through spe- 
cies successional turnover. To quantify the importance of the direct effects and 
the effects mediated by species change in boreal spruce swamp forests, we used 
a dual approach: (i) we measured successional changes in moss communities at 
36 sites (nine undrained, nine drained, 18 rewetted) and (ii) photosynthetic 
properties of the dominant Sphagnum and feather mosses at nine of these sites 
(three undrained, three drained, three rewetted). Drainage and rewetting 
affected moss carbon assimilation mainly through species successional turnover. 
The species differed along a light-adaptation gradient, which separated shade- 
adapted feather mosses from Sphagnum mosses and Sphagnum girgensohnii 
from other Sphagna, and a productivity and moisture gradient, which separated 
Sphagnum riparium and Sphagnum girgensohnii from the less productive S. an- 
gustifolium, S. magellanicum and S. russowii. Undrained and drained sites har- 
bored conservative, low-production species: hummock- Sphagna and feather 
mosses, respectively. Ditch creation and rewetting produced niches for species 
with opportunistic strategies and high carbon assimilation. The direct effects 
also caused higher photosynthetic productivity in ditches and in rewetted sites 
than in undrained and drained main sites. 
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Introduction 

Approximately 15 million hectares of Sphagnum peat- 
lands have been drained to enhance tree growth for for- 
estry, mostly in northern Europe (Joosten and Clarke 
2002). Drainage decreases Sphagnum cover (Laine et al. 
1995; Korpela 2004), leading to cessation of Sphagnum 
biomass accumulation and consequently, to a loss of 



many ecosystem services that Sphagnum provides [e.g., 
filtration of soluble organic matter and nutrients, carbon 
store function, and sustenance of species of conservation 
value (Zak et al. 2011)]. Feather mosses, which have a 
lower ability to accumulate carbon than Sphagnum 
(Turetsky et al. 2010), show an opposite pattern: their 
relative cover increases after drainage (Laine et al. 1995; 
Korpela 2004). 
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Spruce swamp forests are minerotrophic peatlands with 
a dense canopy of trees (Picea abies). Despite liaving liigh 
biodiversity (Romberg et al. 1998) and ecosystem service 
values when undrained, spruce swamp forests have been 
among the peatlands most commonly selected for drain- 
age, due to their high productivity when drained for for- 
estry (Joosten and Clarke 2002). In Finland, where 
drainage for forestry has been most intensive (Paivanen 
and Hanell 2012), the area of undrained spruce swamp 
forests has declined from 2 million hectares in the 1950s 
to 0.8 million hectares (Ilvessalo 1958; Virkkala et al. 
2000). The largest decline has occurred in southern Fin- 
land, where spruce swamp forests are now classified as a 
threatened habitat type (Kaakinen et al. 2008). Restora- 
tion of spruce swamp forests started in Finland in the 
1990s. It is mostly done for nature conservation purposes 
in protected areas: rewetting is accomplished by blocking 
the ditches (Aapala and Tukia 2008). Rewetting practices 
are weD developed by now (Aapala and Simila 2013), but 
ecological restoration success remains to be quantified. 

Restoration success can be defined as when the restored 
site follows a trajectory that leads to pristine-like environ- 
mental conditions, and communities and ecosystem func- 
tions typical of pristine ecosystems (Dobson et al. 1997). 
In Sphagnum peatlands, such as spruce swamp forests, 
this includes restoring the growth of Sphagnum mosses, 
which modify the ecosystem physical conditions and are 
the primary peat-forming plants (van Breemen 1995). 

Restoration by rewetting may affect moss carbon assimi- 
lation directly, and/or indirectly through change in species 
composition. The primary direct effect of rewetting on 
mosses involves a change in moisture: as poikilohydric 
plants that cannot regulate their water uptake and loss, 
mosses are dependent on external moisture. Sphagnum net 
photosynthesis is related to its current moisture content 
that correlates with water table level (Schipperges and Ry- 
din 1998; Strack and Price 2009). Past moisture conditions 
also affect Sphagnum photosynthesis (Schipperges and Ry- 
din 1998). Sphagnum mosses are known to grow well in 
water-saturated conditions (Rochefort et al. 2002), such as 
those that prevail immediately after successful ditch block- 
ing (Aapala and Tukia 2008). Feather mosses grow better 
in wet conditions as well, although their abundance is low 
in wet habitats (Bauer et al. 2007). Feather mosses are nor- 
mally restricted from water-saturated environments 
because of physiological constraints and competitive exclu- 
sion by Sphagnum (MuUigan and Gignac 2001, 2002). 

Rewetting may also affect mosses through change in 
species composition, because species that are specialized to 
different habitats differ from one another in productivity. 
For instance, Sphagnum species of wet microhabitats show 
higher growth rates than species of dry microhabitats 
(Gunnarsson 2005), and feather moss species show lower 
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productivity than Sphagnum (Turetsky et al. 2010). Photo- 
synthetic properties differ between plants that are typical 
to different successional stages: maximum photosynthetic 
capacity, dark respiration and light compensation point 
generally decrease from early- to late-successional species 
(Bazzaz 1979), while the physiological stress experienced 
by the plants increases (Grime 1977). In peatlands, this 
development has been described in the succession from 
early-successional, fast-growing hollow species to drought- 
and irradiance-stressed hummock species (Granath et al. 
2010; Laine et al. 2011b). However, as spruce swamp for- 
ests do not have a true hummock-hollow structure and 
lack the high irradiance that causes stress to mosses in 
open mires (Hajek et al. 2009), the successional pattern is 
likely to be different. 

Ditches constitute a distinct habitat in drained and 
rewetted peatlands. In the drained phase, ditches function 
as a refuge for Sphagnum (Laine et al. 1995). Following 
rewetting. Sphagnum biomass in the blocked ditches can 
help to stabilize site hydrology. Ditches and ditch banks 
differ from the main site in water table levels, received 
irradiance and disturbance regime, which may affect moss 
photosynthetic properties directly or via changes in spe- 
cies composition. 

Understanding the mechanisms of Sphagnum recovery 
is vital for understanding the trajectories that lead to peat- 
land restoration success. Thus far, studies on Sphagnum 
growth traits along primary (Laine et al. 2011b) or sec- 
ondary succession (Granath et al. 2010) have concentrated 
on unforested open mires. In this study, we focus on the 
impacts of drainage and rewetting on moss photosynthesis 
in spruce swamp forests. Measurements on CO2 exchange 
provide information on the photosynthetic efficiency and 
light responses of the mosses, while chlorophyll fluores- 
cence measures levels of plant stress due to water limita- 
tions, light intensity, and/or nutrient supply (Maxwell and 
Johnson 2000). We expect drainage and rewetting to affect 
moss carbon assimilation directly and indirectly by chang- 
ing the moss species composition. Our aim is to quantify 
the importance of the direct effects and the indirect effects 
mediated by successional species change. 

This study focuses on five parameters: (i) the maximum 
rate of light-saturated gross photosynthesis (PMAX) show- 
ing the photosynthetic capacity, (ii) dark respiration (_R), 
(iii) light compensation point of net photosynthesis 
(PPFD^,): a measure of photosynthetic light-use efficiency 
at low light, (iv) actual quantum yield of PSII in high light 
('3'psii) showing the efficiency of the photosynthetic 
machinery, and (v) maximum potential quantum yield of 
PSII (FJFn^): a plant stress indicator. Based on ecological 
knowledge on succession (Grime 1977; Bazzaz 1979) 
presented above, we expect Sphagnum photosynthetic 
capacity (PMAX) to be highest in rewetted sites and in 
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ditches, the early successional habitats; intermediate in 
undrained sites, the mature habitats; and lowest in drained 
sites, the suboptimal habitats. Conversely, we expect plants 
stress levels (measured as decreased FJF^) to be highest in 
drained, intermediate in undrained and lowest in rewetted 
sites and in ditches. We expect respiration to be highest in 
drained and rewetted sites and lowest in undrained sites. 
We expect the light compensation point (PPFDc) to be 
low everywhere but in the ditches, which lack tree cover. 

Methods 

We used a dual approach to quantify the impact of drain- 
age and rewetting on mosses. We measured successional 
changes in moss communities at 36 sites (nine undrained, 
nine drained, 18 rewetted) and photosynthetic properties 
of the dominant Sphagnum and feather mosses at nine of 
these sites (three undrained, three drained, and three rew- 
etted). 



Study sites 

Originally, before drainage, all sites were similar {Vaccini- 
um myrtillus spruce mires, Laine et al. 2012). To enhance 
tree growth for forestry, ditches were constructed between 
1900 and 1980. Drainage had increased tree volume in 
the sites from undrained levels (Table A 1.1 in Appen- 
dix 1, P = 0.45). Rewetting was conducted between 1995 
and 2008 (2001-2003 in the sites sampled for the photo- 
synthesis measurements) by the Finnish state forest 
agency Metsahallitus by blocking the drainage ditches 
with peat or wood (Table 1). 

Sites all have peat depths >80 cm and are located in 
the southern boreal zone, 60— 62°N, 23— 25°E (for a map, 
see Appendix 1, Fig. Al.I), at altitudes of 40-170 m a.s.L. 
Climate is cold and humid with a long-term mean annual 
temperature of 3.5-5.3°C and annual precipitation that 
ranges from 627 to 768 mm depending on location 
(Table Al.l in Appendix 1). The average summer 2011 



Table 1. Moss species sampled by site and drainage state. 



Site 


Drainage state 


Year of rewetting 
(drainage) 


Sampled moss species' 
May-August 




EvLuVK 


Undrained 




P. schr 




S. angu 








S. girg 












S. mage 












S. wulf 






SusiLu 


Undrained 




P. schr 




S. angu 








S. girg 




S. mage 








S. ripa 






EvLuPa 


Undrained 




S. russ 




P. comm 








S. girg 




S. angu 








S. mage 






Ev03ku 


Rewetted (via ditch filling) 


2003 (1949-1980) 


P. schr 


S. girg^ 


H. sple 










S. ripa^ 


S. angu 












S. russ 


EvOBma 


Rewetted (via ditch filling) 


2003 (1949-1980) 


P. schr 


S. ripa^ 


S. angu 








5. girg 


S. russ^ 


S. wulf 


EvOIVR 


Rewetted (via ditch blocking) 


2001 (1949-1980) 


P. schr 


S. ripa^ 


S. angu 








S. girg 




S. russ 












5. wulf 


Lakl<Oj 


Drained 


(1949) 


P. schr 




S. russ 








S. girg 




S. mage^ 








S. mage 






KoniOj 


Drained 


(1965) 


P. schr 




S. angu 








S. girg 




S. russ 








S. mage 






VesiOj 


Drained 


(1908-1913) 


P. schr 


S. ripa^ 


H. sple 








S. girg 







' H. sple = Hylocomium splendens, P. schr = Pleurozlum schreberl, P. comm = Poiytrlchum commune, S. angu = Sphagnum angustifolium, 
5. girg = S. girgensohnii, S. mage = S. magellanicum, S. ripa = S. riparium, S. russ = S. russowii S. wulf = S. wulfianum. 



Additional to the species sampled at all times. 
^Sampled from the ditch. 
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May-August temperature in the sites sampled for the 
photosynthesis measurements was 14.6°C, which is 1.7°C 
warmer than the long-term average (1971-2000). Total 
summer 2011 precipitation was 230 mm, 49 mm less 
than the long-term average. Norway spruce {Picea abies) 
was the dominant overstorey species at all sites; the un- 
derstorey was dominated by Vaccinium spp. dwarf shrubs. 



Moss cover survey 

A vegetation survey was conducted at 36 sites (nine 
undrained, nine drained, 18 rewetted, see Appendix 1) in 
2009. In each site, percent cover of each moss species was 
estimated in a total of 72 sample plots, 30 cm in diame- 
ter, placed in a clustered design. Sphagnum girgensohnii 
and Sphagnum russowii were pooled, because they could 
not be visually identified from each other without exten- 
sive effort. 



Photosynthesis measurements 

Sampling and sample preparation 

We measured photosynthesis of Sphagnum and feather 
mosses monthly during the summer of 2011. The sam- 
pling was designed to account for both the direct and 
indirect effects of drainage and rewetting: drainage state, 
variation related to the presence of the ditch habitat and 
differences between moss species. To eliminate the effect 
of short-term fluctuation in moisture, the measurements 
were conducted on acclimatized, moist moss shoots. 
Dominant moss species in each site (3-4 species, except 
for the first sample date 4-6) were collected from the 
most typical habitat for each species (Table 1). Sphagnum 
girgensohnii (Fig. 1, left) and Pleurozium schreberi (Fig. 1, 




Figure 1. Sphagnum girgensohnii (left) and Pleurozium schreberi 
(right) common mosses in undrained and drained spruce swamp 
forests. Photos: Juicka Laine. 



right), which were common to all sites, were always 
collected regardless of dominance. Mosses were collected 
either from near the ditch ("ditch") or away from the 
ditch ("main site") in the drained and rewetted sites 
(Table I). 

Three replicates per species were collected each mea- 
surement period by cutting the top 5 cm of stems from a 
25 cm^ area. At each moss collection point, peat moisture 
of the top 12 cm was measured using a CS-620 Hydro- 
Sense (Campbell Scientific, Logan, UT) meter. During 
each sampling period, site water table (WT) was mea- 
sured manually from three perforated wells that transect- 
ed the center of each site. In the drained and rewetted 
sites, one well was located in the ditch line, and two wells 
transected the main site. Data on tree volume were 
obtained from tree stand measurements conducted in the 
sites in 2010. Mosses were stored in polyethylene bags to 
maintain moisture. After field collection, they were kept 
in the dark at 5°C for up to 2 days until photosynthesis 
could be conducted. 

CO2 exchange and chlorophyll fluorescence 
measurements 

CO2 exchange and chlorophyll fluorescence were mea- 
sured in the laboratory using a portable gas exchange flu- 
orescence system GFS-3000 (Heinz Walz GmbH, 
Effeltrich, Germany). We used a 4 x 2 cm standard 
chamber that was modified to measure photosynthesis on 
moss shoot segments. The opaque plexiglass cuvettes 
(frames) were 1 cm high and equipped with a mesh bot- 
tom surface to allow free airflow around the sample. A 
uniform layer of Sphagnum capitula (corresponding to 
the top 1 cm) was placed in the cuvette. The number of 
capitula used varied by species and ranged from 5 to 16. 
For feather mosses, the top 2 cm were cut and placed 
lengthwise in the cuvette, with stem numbers ranging 
from 4 to 11. In order to homogenize and minimize 
water content of the shoots arranged in the cuvettes, we 
saturated them with drops of water and then blotted 
them gently from both sides with sheets of pulp until they 
released no more water. We verified that under these 
experimental conditions the shoot water content range 
represented the optimum for CO2 exchange. 

Prior to measurements, the dark- acclimated samples 
were allowed to acclimate in the cuvettes for 20 min 
under a PPFD of 1000 /imol nC^ s^^ and ambient room 
temperature of approx. 22°C. Net photosynthesis (A) was 
measured at decreasing levels of PPFD: 1000, 50, 25, and 
0 /mol m^^ s^' (abbreviated as Aiooo> ^so. ^25> and Aq) 
with artificial light provided by a built-in LED light 
source. A 1000 represents the maximum photosynthetic 
capacity of the mosses, A50 to A25 show net photosyn- 
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thetic rate in shaded conditions and Aq represents respira- 
tion. Samples were allowed to acclimate to each light level 
prior to measurement until A was constant. During the 
measurement period, the chamber temperature was kept 
constant at 20°C, the CO2 concentration of incoming air 
was 400 ppm, air flow was 400 /^mol s^' (9.6 mL s^'), 
and the relative humidity of outgoing air was maintained 
at approximately 90%. Light compensation point of net 
photosynthesis (PPFD^) was defined as the level of 
PPFD where A = 0, calculated from the initial part of the 
A/PPFD curve (from Aq to A50). Because A50 was used for 
deriving PPFDc, it was not used as an independent vari- 
able in the further data analysis. 

Parameters related to photosystem II (PSII) were mea- 
sured to assess the amount of stress experienced by the 
mosses, which reflects acclimation of the mosses to their 
habitats. 

Actual quantum yield of PSII photochemistry ($psii) 
and maximum fluorescence (F^) were measured at the 
end of the 1000 /imol m^^ s^' light level. Samples were 
then dark acclimated for 6-12 h at 5°C. After the dark 
acclimation, chlorophyll fluorescence was measured again, 
and the ratio of variable and maximum fluorescence 
(fv/f^m) calculated. The FJF^-n ratio represents the maxi- 
mum potential quantum yield of PSII. After the measure- 
ment, samples were dried to a constant weight, and A 
was expressed per unit dry mass (mg g^' h^^). 

Data analysis 

In the analysis of the effects of drainage and rewetting, 
we used both classification into drainage states 
(undrained, drained, rewetted), and classification into 
habitats (undrained, drained main site, ditch of drained 
site, rewetted main site, ditch of rewetted site). The lat- 
ter one acknowledges the marked spatial variation 
related to the presence of the ditch. In addition, as pre- 
vious studies have revealed that photosynthetic responses 
of peatland mosses often vary by season (i.e., Gaberscik 
and Martincic 1987), this was taken into account in the 
analysis. 

To quantify the effect of habitat on moss community 
composition in the 36 sites of the vegetation survey, we 
used redundancy analysis (RDA) on centered, nontrans- 
formed moss species data using the program Canoco 5 
(ter Braak and Smilauer 2012). Statistical significance was 
evaluated using Monte Carlo permutation restricted for 
the hierarchical sampling design. 

To quantify differences in water table level between the 
habitats in the nine sites sampled for the measurements, 
we applied a linear mixed-effects model. In the initial 
model, habitat, month and the interaction of these two 
were included as fixed effects. Site was included as a ran- 



dom effect. The interaction was not found significant and 
was eliminated. Differences in water table level between 
the habitats and months were compared post hoc. We 
quantified differences in tree stand volume between the 
drainage states in the nine sites using ANOVA and post 
hoc comparisons. Models were fitted using functions Ime 
and Imer in the lme4 package of R. 

To explore the main trends in the variation of the mea- 
sured photosynthetic response parameters, we used prin- 
cipal component analysis (PCA) on CO2 assimilation rate 
at three levels of PPFD (Aiooo, A25, and Aq), light com- 
pensation point of net photosynthesis (PPFDc), actual 
quantum yield of PSII in high PPFD (<I>psii), and maxi- 
mum potential quantum yield of PSII (FJF^). The varia- 
tion in the parameters was projected in relation to 
habitat, species, peat field moisture, and site water table. 
As patterns without seasonal variation are easier to inter- 
pret, only data from the May measurement period was 
used for the PCA analysis. This month contained the 
greatest number of measured species. 

To quantify the direct effect of drainage and rewetting 
against the effect of moss species on moss photosynthetic 
parameters, we conducted two variation-partitioning 
analyses: one using the drainage state and another one 
using the habitat as a predictor variable. The first analysis 
provides a direct answer to our research question, 
whereas the second one acknowledges the actual habitat 
diversity created by drainage and rewetting. We parti- 
tioned the variation in the measured photosynthetic 
response parameters into three components explained by 
species, month, and drainage state/habitat, testing both 
simple and conditional effects. This was conducted by 
creating a partial RDA for each predictor variable with 
the other predictor variables as covariates. Only the spe- 
cies measured in aU 4 months were included in the RDA. 
CANOCO for Windows 4.5 and 5 (ter Braak and Smila- 
uer 2002, 2012) was used for the PCA and RDA. The 
analyses were conducted on centered and standardized 
photosynthetic parameters as response variables. 

To quantify differences between light compensation 
point of net photosynthesis (PPFD^), actual quantum 
yield of PSII ($psii), and the maximum potential quan- 
tum yield (-Fv/^^m)i we applied linear mixed-effects mod- 
els. In the initial models, species, habitat, month, water 
table, peat field moisture, and sample dry weight were 
included as fixed predictors. Site was included as a ran- 
dom effect. Fixed effects were eliminated from the model 
if not found significant (see Table A2.3 in Appendix 2 for 
the final model results). The differences in PPFDc, ^psn, 
and Fy/Fn^ between the habitats, species, and months were 
compared post hoc as described previously. Models were 
fitted using functions Ime and gls in the nlme package of 
R (Pinheiro and Bates 2000). 
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To determine the effects of habitat and species for pho- 
tosynthetic capacity, light-use efficiency at low light, and 
respiration, we applied a nonlinear mixed-effects model 
with the hyperbolic light saturation curve (i.e., Larcher 
2003), within which parameters were linearly dependent 
on predictor variables: 



PMAXksPPFD^si , 
a + PPFDisi 



(1) 



where the response Ai^i is the observed net photosynthesis 
and the predictor PPFDf^i is the photosynthetic photon 
fltix density for measurement i of sample s on site k. The 
parameters to be estimated are respiration {R^s) photo- 
synthetic capacity i.e. the maximum rate of light-saturated 
gross photosynthesis (PMAXi^) and the maximum quan- 
tum yield of CO2 assimilation (a). The residual (e^;) is 
normally distributed with mean zero and constant vari- 
ance. Parameter a. was assumed to be constant over all 
samples and sites; this restriction was necessary because 
of the low number of measurements per sample (four 
PPFD levels with one observation for each). Variation in 
Rks and PMAXi^ was explained by the fixed predictors 
moss species, habitat, month, water level, peat field mois- 
ture and sample dry mass, and random effects for site 
and sample. Final models for the photosynthesis parame- 
ters in Eq. 1 are defined below (see Table A2.4 in Appen- 
dix 2 for the final model results). AU terms in the 
following models explained the variation in response sig- 
nificantly (approximate F- test, P < 0.05): 



Rks = SPks + MOks + Hfa + rk + ns 



(2) 



PMAXks = SPks + MOks + Hfa + MCks + + flfa (3) 

where SPj^, MO^s, and Hjts are factor-type predictors for 
species (9 levels), month (4 levels), and habitat (5 levels), 
respectively. MQs is dry mass of the sample, which has 
been centralized to have a mean of zero. The last two 
terms in the equations are random effects for the site and 
sample, with bivariate normal distributions (r^, a^Y ~ 
MVN{0, E^) and (r^, a^)' ~ MViV(0, 17^). The random 
effects account for the correlation arising from the nested 
grouping of the data to sites and samples within sites. 
The model was fitted and the tests performed using pack- 
age nlme of the R software (Pinheiro and Bates 2000). 

The differences in PMAX and _R (Eq. 1) between the 
habitats, species, and months were compared post hoc: 
each habitat was compared against undrained, moss spe- 
cies were compared against Sphagnum girgensohnii, and 
months were compared against July. The difference to 
undrained shows how drainage and rewetting have chan- 
ged the photosynthetic parameters from the original natu- 
ral conditions. Of the moss species, S. girgensohnii was 
chosen as the baseline because it is a common, typical 



moss species in undrained spruce swamp forests (Laine 
et al. 2012). July was chosen as the baseline month 
because it is the usual period of peak growth in the study 
region (Riutta et al. 2007; Wilson et al. 2007). Signifi- 
cance limit of P < 0.05 was used in all analyses. 

To test whether the effects of drainage and rewetting 
differ between Sphagnum girgensohnii and Pleurozium 
schreberi - the two species that we sampled in all three 
drainage states - we conducted a separate test on the 
interaction effect of species and habitat on $psii, FJF,„, 
PMAX, and _R. Ditch habitats were excluded from the 
analysis, as P. schreberi did not occur in them. 

To examine photosynthetic strategies of the moss 
species a posteriori, we classified the species in three 
categories after Grime (1977): ruderal, competitive, and 
stress-tolerant, based on their PPFD,- and PMAX. The 
stress-tolerant category was further divided into stress- 
tolerant, shade species; and stress-tolerant, light species. 
Ruderal species were defined to show high PPFD^ and 
PMAX. Competitive species were defined to show low 
PPFD;-, because they are more adapted to the shaded con- 
ditions of the spruce swamp forests, and high PMAX. 
Stress-tolerant species were defined to show low PMAX. 
Stress-tolerant shade species were defined to have low 
PPFDc, stress-tolerant light species high PPFD^. To test 
the classification, we ran the models for the photosynthetic 
parameters using these four groups instead of species. 

Results 

Moss species composition 

Total moss cover was highest in the undrained sites. 
S. girgensohnii coupled with S. russowii was favored by 
undrained conditions, but it was common in all habitats 
(Fig. 2). Pleurozium schreberi and Hylocomium splendens 
were more common in drained and rewetted than in 
undrained sites. High cover of Sphagnum riparium 
and S. squarrosum was typical of ditch habitats, while 
S. magellanicum and S. angustifolium were typical species 
for undrained sites (Fig. 2). 



Photosynthetic properties 

Environmental conditions 

In the sites sampled for the photosynthesis measurements, 
ditches in drained sites had the highest water table, fol- 
lowed by ditches in rewetted sites (Fig. 3A). Rewetted and 
undrained sites showed similar (P = 0.97) water table lev- 
els (Fig. 3A). Water tables were lowest in drained sites, 
but difference to undrained and rewetted sites was not 
significant (Fig. 3A). Differences in water table between 
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Figure 2. (A) Sphagnum and feather moss 
species cover by habitat; (B) redundancy 
analysis (RDA) on the effect of habitat on moss 
community composition. Sphagnum and 
feather moss species with >10% fit shown. 
First axis explains 15% of the data variation, 
P = 0.004. Second axis explains 7% of the 
data variation, P = 0.002. 
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the habitats remained similar over the growing season of 
2011, as indicated by a lack of significant interaction 
effect. Water table levels in May and June were signifi- 
cantly higher than water tables in July and August. Tree 
stand volume was highest in the drained sites (Fig. 3B), 
but not significantly so (0.05 < P-values<0.10). 

Main gradients in the data 

Two strong gradients appear in the photosynthetic 
response data (Fig. 4). The main gradient (PCA Axis 1) is 
related to photosynthetic efficiency at low light/dark res- 
piration (Ao), photosynthesis at the lowest light level 
(A25) and light compensation point of net photosynthesis 
(PPFDc). It separates feather mosses from Sphagnum 



mosses and S. girgensohnii from the remaining Sphagna. 
This gradient explains 44% of the variation in photosyn- 
thetic properties. The second gradient (PCA Axis 2) is 
related to productivity and moisture: photosynthesis at 
high light (Aiooo)> the stress indicator FJF^, water table 
and peat field moisture. It reflects differences in both spe- 
cies and habitats. Along this gradient, increased field 
moisture corresponds to higher productivity and 
decreased stress. The second gradient explains 30% of the 
data variation. 

Species versus habitat influence 

Moss species affected variation in the measured photosyn- 
thetic properties more than drainage state (Table 2a): spe- 



Figure 3. (A) Average water table level (WT) 
relative to moss surface during the summer 
season 201 1 by habitat and (B) tree stand 
volume in the study sites by drainage state; 
Bars indicate SE. Different letters mark 
significant differences in ANOVA; letters in a) 
apply to all months, although marked only for 
July. 
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▲ 
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Ditch of ^ 
drained site S. riparium 
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Figure 4. Principal component analysis (PCA) linking photosynthetic 
response parameters with environmental factors. Data measured 
during May 201 1 . Photosynthetic response parameters: CO2 
assimilation rate at three levels of PPFD Wiooo. ^25. and Aq), light 
compensation point of net photosynthesis (PPFD^), maximum 
quantum yield of PSIl (FJF^), and quantum yield of PSIl (<i>psii). 
Environmental factors: moss species, habitat, site water table (wf), 
and peat field moisture (pfm). Axes 1 (light adaptation) and 2 
(productivity and moisture) explain 44% and 30% of total variation, 
respectively. 

cies explained 31% of the variation, while drainage state 
explained 2%. When the presence of ditches is taken into 
account (Table 2b), combined effect of species and habi- 
tat became important, explaining 7% of the data varia- 
tion: this reflects species differences between ditch 
habitats and the main sites (Table 1, Fig. 2). 

Both species and habitat affected the modeled photo- 
synthetic parameters PPFDc, PMAX and R significantly 
(Tables A2.1 and A2.2 in Appendix 2). Also, ^psn was 
affected by both species and habitat (Table A2.1). Stress 
indicator Fy/F„ was affected by species, but not habi- 
tat: instead, water table level and field moisture were 
significant predictors for Fy/Fm (Table A2.1). Parameter 
estimates and standard errors for the fixed effects, and 
standard deviations and correlations for the random 
effects are presented in Tables A2.3 and A2.4 in Appen- 
dix 2. 



Light adaptation 

Low light compensation point of net photosynthesis 
(PPFD^) is an adaptation to shady environment. PPFD,- 
was lower in undrained than in drained main sites 
(Table 3a). The lowest PPFDc was found in the ditches of 



Table 2. Variation partioning into fractions explained by moss spe- 
cies, month and (a) drainage state (undrained, drained, rewetted) or 
(b) habitat (undrained, drained, ditch of drained site, rewetted, ditch 
of rewetted site). All fractions were significant, P < 0.005. 



Predictor 



% of all % of all 

variation Predictor (combined effects) variation 



(a) 



Species 
Month 

Drainage state 



Total 

(b) 
Species 
Month 
Habitat 



Total 



30.5 
16.9 
1.8 



53.1 

25.4 
17.1 
3.3 



54.5 



Species & Month 2.1 

Species & Drainage state 1.7 

Month & Drainage state <0.1 

Species, Month & Drainage -0.1 
state 



Species & Month 3 
Species & Habitat 6.8 

Month & Habitat -0.1 

Species, Month & -1 
Habitat 



drained sites (Table 3a). Feather mosses Pleurozium schre- 
beri and Hylocomium splendens had the lowest and Sphag- 
num riparium the highest PPFDc (Table 3b). S. wulfianum 
and S. girgensohnii had lower PPFDc than the other Spha- 
gna (Table 3b). 

Productivity 

PMAX and R were higher in ditches and in rewetted sites 
than in undrained sites, but similar across undrained and 
drained main sites (Table 4a). Feather mosses, Pleurozium 
schreberi and Hylocomium splendens, had the lowest and 
Sphagnum riparium the highest PMAX, R and net produc- 
tivity (Table 4b). S. russowii and S. magellanicum had 
lower PMAX than S. girgensohnii but similar R, which 
resulted in lower net maximum productivity than that of 
S. girgensohnii (Table 4b). 

Physiological efficiency and stress 

Ditches of drained sites had lower <I>psii than the other 
habitats (Table 3a). For H. splendens, P. commune, and 
P. schreberi, ^psn was 50% higher than for Sphagnum 
mosses (Table 3b). F^/F^ responded to water level and 
field moisture, not to habitat (Table A2.1). F^/Fm was low- 
est, i.e. stress was highest, for S. riparium and S. magellan- 
icum; Polytrichum commune showed the highest Fy/F^ 
(Table 3b). 

Seasonality in photosynthetic properties 

PMAiX showed no change across the season (Table 4c), 
but R (Table 4c) and PPFD^ (Table 3c) were higher in 
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Table 3. Differences in light compensation point {PPFD^), maximum quantum yield of PSIl (FJFm), and quantum yield of PSIl ($ps/;) between (a) hab- 
itat, (b) species (in the order of increasing PPCDJ, and (c) month. Posf hoc contrast results from the linear-mixed-effects models. Undrained state, 
Sphagnum girgensohnii and July are the baselines, two of which are kept constant while the predictor variable in question changes. P-values 
indicate significant differences from undrained, S. girgensohnii, and July, respectively. Average ± SE, n = number of measured moss samples. 





n 


PPFD^ (/imol m"^ s"'; 


1 P-value 


$PS|| 




P-value* 




(a) Habitat S. girg, July 
















Undrained 


49 


14 ± 1 


- 


0.09 ± 0.004 


- 




Drained, main site 


35 


16.5 ± 1.19 


0.035 


0.09 ± 0.003 


0.596 




Rewetted, main site 


29 


16 ± 1.21 


0.096 


0.1 


± 0.004 


0.145 




Drained, ditch 


6 


9.9 ± 1.95 


0.034 


0.06 ± 0.005 


<0.001 




Rewetted, ditch 


19 


14.7 ± 1.52 


0.655 


0.09 ± 0.003 


0.838 






n 


PPFD, inmol s"') 


P-value* 


$PSII 


P-value* 




P-value 


(b) Species Undrained, July 
















Hyiocomiunn spiendens 


2 


8.83 ± 2.24 


0.022 


0.16 ± 0.02 


0.001 


0.76 ± 0.012 


0.7 


Polytrichum commune 


1 


11.3 ± 3 


0.361 


0.16 ± 0.028 


0.023 


0.82 ± 0.016 


<0.001 


S. wulfianum 


6 


13.1 ± 2.2 


0.684 


0.12 ± 0.007 


<0.001 


0.77 ± 0.007 


0.111 


Pieurozium schreberi 


36 


14 ± u.y/ 


0.998 


0.17 ± 0.005 


<0.001 


U./D ± U.UU4 


0.001 


S. girgensohnii 


36 


14 ± 1 




0.09 ± 0.004 




0.76 ± 0.007 




S. magelianicum 


18 


18.3 ± 0.96 


<0.001 


0.09 ± 0.003 


0.275 


0.73 ± 0.005 


<0.001 


S. angustifolium 


7 


19 ± 1.88 


0.009 


0.11 ± 0.007 


0.01 


0.75 ± 0.007 


0.109 


S. russowii 


12 


20.5 ± 1.27 


<0.001 


0.1 ± 0.004 


0.006 


0.76 ± 0.005 


0.481 


S. riparium 


20 


22.4 ± 1.24 


<0.001 


0.1 ± 0.003 


0.353 


0.72 ± 0.006 


<0.001 


(c) Month S. girg, Undrained 
















May 


48 


24.1 ± 0.79 


<0.001 


0.12 ± 0.003 


<0.001 


0.738 ± 0.004 


<0.001 


June 


30 


15.4 ± 0.65 


0.04 


0.1 ± 0.002 


0.194 


0.775 ± 0.004 


<0.001 


July 


30 


14 ± 1 




0.09 ± 0.004 




0.76 ± 0.007 




August 


30 


14.8 ± 0.64 


0.218 


0.11 ± 0.003 


<0.001 


0.802 ± 0.004 


<0.001 



Bold font indicates relationship is significant. 



May and June than in July and August. Plant stress, as 
indicated by low FJF^, was highest in May and lowest in 
August (Table 3c). 

Habitat effect by species 

Sphagnum girgensohnii and Pieurozium schreberi did not 
differ in their responses to drainage state in most photo- 
synthetic parameters. Only the plant stress indicator 
Fv^Pm, showed a larger difference for the drained state to 
undrained and rewetted states for S. girgensohnii than 
for Pieurozium schreberi (P = 0.019). The Fy/F„, values 
(lower values for higher stress) for P. schreberi in 
undrained, drained and rewetted conditions were 0.74, 
0.74 and 0.75, respectively; for S. girgensohnii 0.76, 0.74 
and 0.77. 



Moss strategies 

The four groups: stress-tolerant (shade), ruderal, competi- 
tive and stress-tolerant (light) (Table 5) functioned as sig- 
nificant predictors in the models for PPFD^ and PMAX 
when used as substitutes for species (see Fig. 5 for the 
estimates). 



Discussion 

Spruce swamp forest 
mosses 



favorable habitat for 



Minerotrophic, shaded mire environments provide favor- 
able growing conditions for Sphagnum mosses (Clymo 
1973; Hajek et al. 2009). These conditions are found in 
rewetted and undrained spruce swamp forests throughout 
the growing season. Values for PMAX were higher than 
those reported for Sphagnum and feather mosses in om- 
brotrophic bogs (Granath et al. 2009; Laine et al. 2011b), 
forested permafrost peatlands (Skre and Oechel 1981), a 
rich fen (Granath et al. 2009) and oligotrophic fens (La- 
ine et al. 2011b). Photosynthetic rates rose gradually from 
spring to mid-season (July), in contrast to ombrotrophic 
bogs, where moss growth tends to be greatest in the 
spring and late summer or autumn (Silvola and Heikki- 
nen 1979; Lindholm 1990; Laine et al. 2011b). A similar 
gradual rise and mid-season peak in photosynthesis has 
been measured in a black spruce permafrost peatland in 
interior Alaska (Skre and Oechel 1981). 

Although water table levels were progressively lower 
toward late summer, values of F^/F^ revealed no drought 
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Table 4. Differences in maximum pinotosynthetic rate (PMAX) and dark respiration (R) between; (a) habitat, (b) species (in the order of increasing 
PMAX) and (c) month. Posf hoc contrast results from the hyperbolic light saturation model (Eq. 1). Undrained state, Sphagnum girgensohnii and 
July are the baselines, two of which are kept constant while the predictor in question changes. P-values indicate significant differences from 
undrained, 5. girgensohnii and July, respectively. Average ± SE, n = number of measured moss samples. 

PMAX+ R 

n PM/\X (mg g"' h"') P-value* R (mg g"H"') P-value* (mg g"H"^) 



(a) Habitat S. girg, July 



Undrained 


49 


6.73 ± 0.31 




-0.831 


± 


0.075 




5.9 


Drained, main site 


35 


6.78 ± 0.42 


0.902 


-0.950 


± 


0.051 


0.019 


5.8 


Rewetted, main site 


6 


7.40 ± 0.42 


0.113 


-1.068 


± 


0.092 


0.01 


6.3 


Drained, ditch 


19 


7.71 ± O.BB 


0.076 


-0.794 


± 


0.131 


0.778 


6.9 


Rewetted, ditch 


29 


7.89 ± 0.47 


0.014 


-1.120 


± 


0.107 


0.007 


6.8 


(b) Species Undrained, July 


















Pleurozium schreberi 


36 


2.68 ± 0.24 


<0.001 


-0.132 


± 


0.051 


<0.001 


2.5 


l-lylocomium splendens 


2 


3.21 ± 0.61 


<0.001 


-0.010 


± 


0.159 


<0.001 


3.2 


S. russowii 


12 


5.46 ± 0.29 


<0.001 


-0.864 


± 


0.075 


0.662 


4.6 


S. magellanicunn 


18 


5.55 ± 0.24 


<0.001 


-0.837 


± 


0.064 


0.923 


4.7 


S. angustifolium 


7 


6.54 ± 0.34 


0.569 


-0.996 


± 


0.091 


0.071 


5.5 


5. girgensohnii 


36 


6.73 ± 0.31 




-0.831 


± 


0.075 




5.9 


Poiytrichum commune 


1 


7.84 ± 0.83 


0.184 


-0.923 


± 


0.22 


0.677 


6.9 


S. wulfianum 


6 


8.01 ± 0.45 


0.004 


-0.982 


± 


0.098 


0.125 


7 


S. riparium 


20 


8.7 ± 0.3 


<0.001 


-1.645 


± 


0.078 


<0.001 


7.1 



PMAX+ R 

n R (mg g"' h"') P-value* (mg g"' h"^) 



(c) Month S. girg, Undrained 



May 


48 


-1.369 ± 0.041 


<0.001 5.4 


June 


30 


-0.934 ± 0.043 


0.016 5.8 


July 


30 


-0.831 ± 0.075 


5.9 


August 


30 


-0.886 ± 0.043 


0.194 5.8 


Bold font indicates relationship is significant. 




Table 5. Species classified by their 


light adaptation, productivity and strategy, based on the photosynthetic response parameters PPFDc (light 


adaptation) and PMAX (productivity). 










Light adaptation 




Species 




(shade/light) 


Productivity (+/-) Strategy (after Grime 1977) 


Pleurozium schreberi 




Shade 


- Stress-tolerant (shade) 


l-lylocomium splendens 




Shade 


— Stress-tolerant (shade) 


Poiytrichum commune 




Shade 


-1- Competitive 


Sphagnum girgensohnii 




Shade 


-1- Competitive 


S. wulfianum 




Shade 


-1- Competitive 


S. riparium 




Light 


-1- Ruderal 


S. magellanicum 




Light 


— Stress-tolerant (light) 


S. russowii 




Light 


— Stress-tolerant (light) 


S. angustifolium 




Light 


- Stress-tolerant (light) 



stress to photosystem II values. On the contrary, FJFm 
increased toward August. The FJF^^^ values were higher 
than those measured in bryophytes from other natural 
conditions (Hajek et al. 2009; Laine et al. 2011b; Zona 
et al. 2011), close to values measured in unstressed vascu- 
lar plants and mosses (Proctor 2010), indicating low levels 
of light-induced stress. The significant differences we 
found in F^/F^ between habitats, species and species 



responses to drainage states were too small to be ecologi- 
cally relevant. 

Direct habitat effects 

Although drainage for forestry deteriorates the conditions 
for Sphagnum mosses, as the decreased Sphagnum cover 
indicates, some microsite areas in the drained sites remain 
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Figure 5. Results from statistical testing on the groups in Table 5, 
when used in the mixed-effect models as substitutes for species: 
differences in light compensation point (PPFD^) and maximum 
photosynthetic rate (PMAX). st(s) = stress-tolerant (shade), 
CO = competitive, ru = ruderal, st(l) = stress-tolerant (light). Different 
letters mark significant differences. 



suitable for Sphagnum: photosynthetic capacity and net 
production in the mosses of these remnant patches did 
not differ from undrained conditions. Sphagnum mosses 
were slightly more productive in rewetted than in 
undrained and drained conditions, but the largest differ- 
ences occurred between the ditches and other habitats. 
This direct habitat effect worked to the same direction as 
the effect mediated by a change in species composition. 

Species responses reflect their growth 
strategies 

Species differences explained the most variation in the 
photosynthetic response parameters PPFDc, PMAX, and R. 
Photosynthetic responses of species reflect their growth 
strategies. Similarly, to moss species along a primary suc- 
cession chronosequence of mires (Laine et al. 2011b), 
responses of the spruce swamp forests moss species reflect 
environmental gradients in light and moisture. The moss 
species can be classified in the three groups defined by 
Grime (1977) as (i) ruderal species that show high produc- 
tion and occupy recently disturbed areas, (ii) competitive 
species that show high production and occur in more sta- 
ble conditions, and (iii) stress-tolerant species that show 
lower production but are more adapted to stress or 
resource scarcity (Table 5). 

Sphagnum riparium is most commonly found at the 
surface water level (Gignac et al. 1991) and is frequently 
a pioneer species in peatlands that experience a rise in 
water table level (Zoltai 1993). S. riparium displayed 
characteristics of ruderal vascular plants (Grime 1977; 
Bazzaz 1979), with high net productivity, PMAX, R, and 
PPFD^. 

Sphagnum girgensohnii is the dominant moss species in 
rewetted sites and appears to be competively superior. It 
had the highest net photosynthesis of all species in rewett- 
ed sites outside the ditch line. Previous research indicates 



S. girgensohnii to be an opportunist species in new habi- 
tats and a key driver of paludification of boreal maritime 
forests in North America (Noble et al. 1984; Asada et al. 
2004). Similar factors, disturbances to the forest floor 
together with increased water table, contributed to the 
increased dominance of S. girgensohnii in those forests 
and in our rewetted sites. Values of FJFm were always 
high, except for a slight decline in drained sites: the low 
stress level indicates fairly large ecological amplitude for 
this species. Sphagnum girgensohnii also differed from the 
remaining Sphagnum mosses by its lower light compensa- 
tion point, which indicates suitability to the shaded habi- 
tat of spruce swamp forests. 

Feather mosses P. schreberi and H. splendens had low 
carbon assimilation and dark respiration rates and low 
light compensation points. They could be classified as 
stress-tolerant species, as they are adapted to shaded, dry 
forest conditions. Another group of stress-tolerants is the 
hummock- Sp/ia^Mfl: S. magellanicum, S. russowii, and 
S. angustifolium. They are not specifically adapted to the 
shaded conditions of spruce swamp forests but tolerate 
drought by forming tight cushions (Clymo 1973). 

Change in ecosystem photosynthetic traits 
through succession 

Species turnover along the sequence of changed conditions 
- drainage and ditch creation, rewetting, and development 
of pristine-like conditions - affects photosynthetic proper- 
ties of the spruce swamp forest ecosystem. Both undrained 
and drained spruce swamp forests can be compared with 
the late-successional stage of forested vascular plant com- 
munities, where succession is associated with decreased 
availability of resources (Grime 1977). Hummock-Sphagna 
(S. magellanicum, S. russowii, and S. angustifolium) are 
typical species of the undrained late-successional stage, 
while feather mosses (Pleurozium schreberi and Hylocomi- 
um splendens) are typical of the drained late-successional 
stage. Ditch creation and rewetting are disturbances that 
create niches for species with opportunistic strategies. In 
our study, main PCA gradients separated the drier, more 
stabilized undrained and drained sites from the wet and 
disturbed rewetted sites and ditch habitats. The three spe- 
cies strategies, as defined by Grime (1977), can be placed 
along the successional gradient: stress-tolerant P. schreberi, 
S. magellanicum, S. russowii, and S. angustifolium at the 
late-successional stages, ruderal S. riparium occupying 
recently disturbed areas and competitive S. girgensohnii 
during mid-succession. 

Ditches of drained sites offered a suitable refuge for 
Sphagnum species to persist. In the rewetted sites, ditches 
are habitats of highly productive Sphagnum cover, pri- 
marily S. riparium. Over time, the high rate of production 
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of the ruderal S. riparium will accelerate terrestrialization 
of the ditch line, which will lower the relative water table 
and create suitable microhabitat for other species. 

Implications 

Functional trait analysis is a useful method for assessing 
the outcome of ecological restoration (Hedberg et al. 
2013), but the established traits, developed for vascular 
plants, do not reflect bryophyte ecology or performance 
(Rice et al. 2008). Photosynthetic properties of mosses are 
directly linked to their evolutionary strategies. If they are 
species-specific, as we here show, they can be used as 
traits in functional trait analysis when coupled with plant 
cover estimations. Light compensation point for net pho- 
tosynthesis {PPFDc), and maximum photosynthesis 
{PMAX) appeared useful in understanding the functional 
variation in spruce swamp forest mosses. 

Peatland restoration monitoring commonly depends 
on comparing restored to pristine sites, which implies 
straightforward directional change. This can be justified 
in ombrotrophic bogs, where vegetation changes after 
drainage and rewetting can be small, because few species 
are able to live in such acid and nutrient-poor condi- 
tions (Laine et al. 2011a). In minerotrophic peatlands, 
development after rewetting involves more species turn- 
over along the successional trajectory (Haapalehto et al. 
2011; Hedberg et al. 2012). Although the species and 
trait composition of the rewetted sites differs from 
undrained systems, especially in the blocked ditches, the 
ruderal and competitive species are likely to contribute 
to the rapid biomass production and peat formation 
during the initial stages after rewetting. Later, the ruder- 
al species are likely to become outcompeted by other 
species of Sphagnum. Measurement-based information 
on species functional traits along successional trajectories 
enables restoration monitoring to identify different stages 
of restoration succession. 
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Appendix 1 




Figure A1.1. Location of the spruce swamp sites used for the 
vegetation survey and (circled) for the photosynthesis measurements. 



Table A1.1. Spruce swamp sites used for the vegetation survey and (colored) for the photosynthesis measurements: time of drainage and rewett- 
ing, location as coordinates (EUREF, ~WGS84), mean annual temperature and annual precipitation in the nearest weather station (1971-2000, 
Finnish Meteorological Institute), average water table depth (WT) as centimeters below moss surface (from manual measurements in July-August 
2009, May-June 2010, May 201 1, September 201 1, and May 2012), tree stand volume and Sphagnum cover. 













Mean annual 


Annual 




Tree stand 






Drainage 


Year of 


Year of 




temperature 


precipitation 


Average 


volume 


Sphagnum 


Code 


state 


drainage 


rewetting 


W-coordinate 


(°C) 


(mm) 


WT (cm) 


(m^) 


cover (%) 


SiLuE 


Undrained 






6,686,952 


5.3 


682 


-13 


128 


62 


SiLuW 


Undrained 






6,686,925 


5.3 


682 


-18 


367 


13 


TeLu 


Undrained 






6,683,434 


5.7 


768 


-28 


216 


26 


RuOjSP 


Drained 


1932 




6,692,132 


5.3 


682 


-33 


553 


6 


RuOjSu 


Drained 


1926 




6,693,212 


5.3 


682 


-72 


329 


0 


TeOj 


Drained 


? 




6,684,069 


5.7 


768 


-39 


169 


23 


Nu97 


Rewetted 


1960s 


1997 


6,689,606 


4.6 


647 


-16 


244 


32 


NuOlhi 


Rewetted 


1960s 


2001 


6,687,472 


4.6 


647 


-19 


356 


41 


NuOIW 


Rewetted 


1960s 


2001 


6,687,779 


4.6 


647 


-6 


133 


27 


NuOSku 
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1960s 


2005 


6,689,992 


4.6 


647 


-30 


218 
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NuOSma 
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1960s 


2005 


6,689,683 


4.6 


647 


-5 


237 


11 


NuOSPo 
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1960s 
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6,686,957 
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647 


-13 


319 
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Undrained 






6,799,071 
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645 
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35 


EvLuPa 


Undrained! 
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f -28 
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56 g 


EvLuVK 
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6,791,370 
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47 ■ 
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-7 
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Li98 
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6,733,287 
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-4 
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state 
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rewetting 
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WT (cm) 
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EvOIVR 


Rewetted 


1949-1980 


2001 


6,790,027 
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1965 1 




6,854,362 


3.5 ^^^^1 
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Appendix 2: Statistical models 



Table A2.1. ANOVA results of the linear mixed-effects models for the differences in light compensation point (PPFD^), actual quantum yield of 
PSIl in high light ($ps;;) and maximum potential quantum yield of PSIl (FJF^). WT denotes water table level, VWC peat volumetric water content. 



PPFD, $ps,i FJF, 



Source 


num. df 


den. df 


F-value 


P- value 


num. df 


den. df 


f-value 


P-value 


num. df den. df 


f-value 


P-value 


Intercept 


1 


387 


1331 


<0.001 


1 


387 


50157 


<0.001 


1 391 


13618 


<0.001 


Species 


8 


387 


69 


<0.001 


8 


387 


13 


<0.001 


8 391 


15 


<0.001 


Month 


3 


387 


70 


<0.001 


3 


387 


70 


<0.001 


3 391 


80 


<0.001 


Habitat 


4 


387 


6 


<0.001 


4 


387 


6 


<0.001 








WT 










1 


387 


8 


0.005 


1 391 


21 


<0.001 


VWC 










1 


387 


6 


0.017 


1 391 


6 


0.012 


Dry mass 


1 


387 


46 


<0.001 

















TableA2.2. ANOVA results of the hyperbolic light saturation model (Eq. 1) for the differences in maximum photosynthesis (PMAX) and dark 
respiration (fi). 





num. df 


den. df 




F-value 


P-value num. df 


den. df 


F-value 


P-value 


a 

Constant 


1 


1206 




3223 


<0.0001 








PMAX 










R 








Intercept 
Species 
Habitat 
Dry mass 


1 1206 
8 1206 
4 1206 
1 1206 




464 
49 
2 

162 


<0.0001 
<0.0001 
0.049 
<0.0001 


Intercept 1 
Species 8 
Month 3 
Habitat 4 


1206 
1206 
1206 
1206 


123 
56 
73 
3 


<0.0001 
<0.0001 
<0.0001 
0.022 
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Table A2.3. Parameter estimates of the linear mixed-effects models for the differences in light compensation point (PPFDJ, actual quantum yield 
of PSIl in high light (<I>P5,/) and maximum potential quantum yield of PSIl (FJF^). 

PPFD^ (/imol s"') Sps,, FJF„, 

Coeff. SE P-value Coeff. SE P-value Coeff. SE P-value 



Fixed part 



Constant (S. g. Jul. Prist.) 


14.01 


0.99 


0.000 


0.093 


0.004 


0.000 


0.760 


0.007 


0.000 


Hylocomium splendens 


-5.18 


2.24 


0.022 


0.063 


0.020 


0.001 


-0.005 


0.012 


0.700 


Polytrichum commune 


-2.74 


2.99 


0.361 


0.065 


0.028 


0.023 


0.063 


0.016 


0.000 


S. wulfianum 


-0.90 


2.21 


0.684 


0.026 


0.007 


0.000 


0.012 


0.007 


0.111 


Pleurozium schreberi 


-0.00 


0.97 


0.998 


0.080 


0.005 


0.000 


-0.014 


0.004 


0.001 


S. magellanicum 


4.31 


0.96 


0.000 


-0.003 


0.003 


0.275 


-0.029 


0.005 


0.000 


S. angustifolium 


4.95 


1.88 


0.009 


0.018 


0.007 


0.010 


-0.011 


0.007 


0.109 


S. russowii 


6.45 


1.28 


0.000 


0.011 


0.004 


0.006 


-0.003 


0.005 


0.481 


S. riparium 


8.35 


1.24 


0.000 


0.003 


0.003 


0.353 


-0.037 


0.006 


0.000 


Drained, ditch 


-4.16 


1.95 


0.034 


-0.029 


0.005 


0.000 








Rewetted, ditch 


0.68 


1.52 


0.655 


0.001 


0.003 


0.838 








Rewetted, main site 


2.03 


1.22 


0.096 


0.005 


0.004 


0.145 








Drained, main site 


2.53 


1.19 


0.035 


-0.001 


0.003 


0.596 








May 


10.1 


0.79 


0.000 


0.029 


0.003 


0.000 


-0.022 


0.004 


0.000 


June 


1.34 


0.65 


0.040 


0.003 


0.002 


0.194 


0.015 


0.004 


0.000 


August 


0.80 


0.65 


0.218 


0.016 


0.003 


0.000 


0.042 


0.004 


0.000 


Water table 








0.000 


0.000 


0.000 


0.00049 


0.00011 


0.000 


VWC 














0.00018 


0.00007 


0.012 


Dry mass 


0.11 


0.016 


0.000 














(andom part 




















sd (constant|site) 


1.24 












0.0107 






Residual standard error 


0.94* 


PPFDc'\0.704 




1.04*$psii'\ 


1.78 




0.0264 







Table A2.4. Parameter estimates of the hyperbolic light saturation 
model (Eq.1) based on nonlinear mixed-effects model fit. 



Fixed part Coeff. SE P-value 



F 

Constant 


73.515 


1.294 


0.000 










Constant (S. girgensohnii, July, 


6.733 


0.313 


0.000 


Undrained) 








Pleurozium schreberi 


-4.052 


0.235 


0.000 


Poiytrichum commune 


1.104 


0.831 


0.184 


S. angustifolium 


-0.195 


0.342 


0.569 


S. magellanicum 


-1.186 


0.244 


0.000 


S. riparium 


1.967 


0.302 


0.000 


S. russowii 


-1.275 


0.286 


0.000 


S. wulfianum 


1.273 


0.447 


0.004 


Hylocomium splendens 


-3.526 


0.610 


0.000 


Drained, ditch 


0.979 


0.550 


0.076 


Rewetted, main site 


0.669 


0.422 


0.113 


Drained, main site 


0.051 


0.416 


0.902 


Rewetted, ditch 


1.161 


0.473 


0.014 


Dry mass 


-0.037 


0.00288 


0.000 


; (mg g"' h"') 








Constant (S. girgensohnii. July, 


-0.831 


0.075 


0.000 


Undrained) 








Pleurozium schreberi 


0.699 


0.051 


0.000 


Polytrichum commune 


-0.092 


0.220 


0.677 


S. angustifolium 


-0.165 


0.091 


0.071 



Table A2.4. Continued. 



Fixed part 


Coeff. 


SE 


P-value 


S. magellanicum 


-0.006 


0.064 


0.923 


S. riparium 


-0.814 


0.078 


0.000 


S. russowii 


-0.033 


0.075 


0.662 


S. wulfianum 


-0.151 


0.098 


0.125 


Hylocomium splendens 


0.821 


0.159 


0.000 


June 


-0.103 


0.043 


0.016 


August 


-0.055 


0.043 


0.194 


May 


-0.538 


0.041 


0.000 


Drained, ditch 


0.037 


0.131 


0.778 


Rewetted, main site 


-0.237 


0.092 


0.010 


Drained, main site 


-0.119 


0.090 


0.019 


Rewetted, ditch 


-0.289 


0.107 


0.007 


Random part 


Site 


Sample in site 


sd(PMAX) 


0.451 




1.19 


sd(/?) 


0.092 


0.11 


con{PMAX*R) 


-0.843 


— 1 


0.998 


Residual standard error 


0.444 
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